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Dynamical friction force exerted on spherical bodies 
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ABSTRACT 

We present a rigorous calculation of the dynamical friction force exerted on a spherical 
massive perturber moving through an infinite homogenous system of field stars. By 
calculating the shape and mass of the polarization cloud induced by the perturber 
in the background system, which decelerates the motion of the perturber, we recover 
Chandrasekhar's drag force law with a modified Coulomb logarithm. As concrete ex- 
amples we calculate the drag force exerted on a Plummer sphere or a sphere with the 
density distribution of a Hernquist profile. It is shown that the shape of the perturber 
affects only the exact form of the Coulomb logarithm. The latter converges on small 
scales, because encounters of the test and field stars with impact parameters less than 
the size of the massive perturber become inefficient. We confirm this way earlier results 
based on the impulse approximation of small angle scatterings. 
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1 INTRODUCTION 

The determination of the dynamical friction force exerted by 
a gravitating system of background stars on a test star mov- 
ing through the system is one of th e classical problems o f 
stellar dynamics. In his seminal pap er [Chandrasekharl (| 19431 ) 
envisaged the scenario of a sequence of consecutive gravita- 
tional two-body encounters of t est and field stars in order 
to calculate the drag force (cf. Illenonl Il973l for a modern 
pre sentation ) . In a n alternative approach Marochnikl (Il968l ) 
and lKalnaisI (|l972h determined the collective response of the 
background medium to a moving massive perturber which 
manifests itself as a polarization cloud decelerating the per- 
turber. The dynamical friction force law derived in this way 
is precisely identical to Chandrasekhar's force law, if the 
mass of an individual field star can be neglected against the 
mass of the perturber. The latter is to be expected whenever 
the background can be treated as a gravitating continuum. 

The applications of Chandrasekhar's formula covered a 
wide field, ranging from calculations of the dynamical fric- 
tion and diffusion coefficients of the Fokker-Planck equa- 
tion, which is used to describe the internal dynamics of star 
clusters, to determinations of the sinking rates of satellite 
galaxies which are accreted by massive galaxies. The liter- 
ature on these subjects has become so extensive over the 
years that we do not attempt any detailed review here, be- 
cause this would be beyond the scope of this paper. In the 
context of the sinking satellite problem most of the efforts 
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concentrated on adapting the original concept of an infinite 
homogeneous background of field stars on isotropic straight- 
line orbits to geometries and kinematics appropriate for the 
discs and haloes of spiral galaxies. Moreover, the finite size 
of the perturbing body was taken into account. IWhitej (|l976l ) 
followed Chandrasekhar's concept of a sequence of two-body 
encounters of field and test stars and considered a sequence 
of deflections of field stars by a moving spherical body. Tech- 
nically the deflections were calculated using the impulse ap- 
proximation of small angle scatterrings. The principal result 
was a modification of the Coulomb logarithm so that it does 
not diverge anymore at small scales, because gravitational 
encounters at impact parameters smaller than the size of the 
perturbing body become ineffective. 



In this paper we present a rigorous calculation of the 
dynamical friction force exerted by an infinite homoge- 
nous background on a sph erical mass i ve bo dy u sing the 
wave- mechanical method of iMarochnikl l|l968l ) and iKalnaid 
(|l972l ). This approach is topically closely related to stud- 
ies of the exch ange of angular and linear momentum in 
stellar systems (iLvnden-Bell fc Kalnai jfl973 , |Pekkerlll974l 
Tremaine fc We inberg 1984, Fuchs 2004) or in plasmas (|Stbd 
19921) . It is shown in section 4 below that the shape of the 
perturber affects only the exact form of the Coulomb log- 
arithm. As concrete examples we calculate the drag force 
exerted on a Plumme r sphere and on a sphere with the den- 
sity distribution of a iHernquisd l|l990l ) profile, respectively, 
and compare this with the drag force exerted on a point 
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2 INDUCED POLARIZATION CLOUDS 

We assume an infinite homogenous distribution of field stars 
on isotropic straight-line orbits. The response of the system 
of background stars to the perturbation due to a massive 
perturber is determined by solving the linearized Boltzmann 
equation 



dfx 
dt 



E 



9fi <9£i dfo = Q 
9xi <9xi dvi 



(1) 



where <3?i denotes the gravitational potential of the per- 
turber and /i is the induced perturbation of the distribu- 
tion function of the field stars in phase space. The unper- 
turbed distribution function is described by fo- The solution 
of the Boltzmann equation is greatly facilitated by consid- 
ering Fourier transforms of the perturbations of the distri- 
bution function and the potential, respectively, 



Um; <*Vk expi[ajt + k • x] . 



(2) 



where u and k denote the frequency and wave vector of the 
Fourier components. Without loss of generality the spatial 
coordinates xi and the corresponding velocity components 
Vi can be oriented with one axis parallel to the direction of 
the wave vector. The Boltzmann equation ([TJ takes then the 
form 



ov 







(3) 



with k = |k| and v denoting the velocity component parallel 
to k. Equation (|3]) has been integrated over the two velocity 
components perpendicular to k. In the following we assume 
for the field stars always a Gaussian velocity distribution 
function, or specifically in eq. ([3]) 

dfo 
dv 



v n b 
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(4) 



where rib denotes the spatial density of the field stars. We 
find then the solution of the Fourier transformed Boltzmann 
equation 

. kv 
7w,k = — - 



fib 



— - ^$^,k. (5) 

LO + kv ^2lTO i 

Integrating eq. (JSJ over the u-velocity leads to the den- 
sity distribution of the induced polarization cloud. This has 
been calculated here without takin g into accou nt the self- 
gravity of the background medium. Fuchs ( 2004) has shown 
that in linear approximation the effects of self-gravity can 
be described by another linearized Boltzmann equation of 
the form of eq. (JTJ where $1 denotes then the gravitational 
potential of the density perturbation of the background 
medium. In a self-gravitating system the density perturba- 
tions are the sources of the potential perturbations so that 
the density - potential pair has to fulfill the Poisson equa- 
tion. The solution of the combined Boltzmann and Poisson 
equations describes simply the Jeans collapse of the back- 
ground medium on scales larger than the Jeans length. In 
real stellar systems or dark haloes their Jeans length will 
be always larger than the size of the system, because other- 
wise the system would have collapsed to smaller sizes. Since 
the polarization cloud is contained within the system, self- 
gravity is not important for its dynamics. 

Of course real self-gravitating systems are not ho- 
mogeneous, but their density falls off radially. Moreover, 



their phase space distribution is more complicated than an 
isothermal Gaussian velocity distribution. In such systems a 
moving perturber can incite large scale perturbations which 
contribute to the dynamical friction of the pertu rber as well 
jTremaine fc Weinberg! [l984 IColpi et alJll999h . Obviously 
these effects cannot be described by the s implified model 
adopted here. However, IColpi et al.l fl999) among others 
have shown, that the drag by the localized polarization cloud 
on the moving perturber, which we treat here, is the princi- 
pal effect. 



3 POTENTIALS OF THE PERTURBING 
BODIES 

For reference reasons we include in our analysis also the 
potential of a point mass, which moves with the velocity vo 
along the y-axis, 

Gm 



*i =- 



v / x 2 + (y _ Vot )2 + , 



(6) 



Its Fourier- Transfor m can be calculated using formulae 
(3.754) and (6.561) of lGradshtevn fc Rhvzikl <|2000f > as 



(7) 



(8) 
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Next, the potential ((6} is generalized to 
Gm 

$1 = ; = , 

^Jr% + x 2 + (y - vo*) 2 + z 2 

which corresponds to an exten ded body with the mass di s- 
tribution of a Plummer sphere (jBinnev fc Tremainelll987l ). 

p=i£rfi + 4y f - (9) 

f r o V r oJ 

The Fourier transform of a moving Plummer sphere is 
given by 



Gm To -ik v vni 

$k = ~2^¥^ l(fcro)e 



(10) 



where K\ denotes the modified Bessel function of the second 
kind. As third example we consid er a perturber w hich has 
the mass density distribution of a lHernquistl (|l990l ) profile, 
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2n r(ro + r) 3 
Its gravitational potential is given by 
Gm 

$1 



(11) 



(12) 



ro + r 

The Fourier transform of a m oving Hernquist sphere can b e 
calculated using eq. (3.722) of lGradshtevn fc Rhvzikl jjOQtjfl 
leading to 

r l + kro cos(fcro)si(fcro) 



4> 



2tt 2 k 2 ' 
— kro sin(fero)ci(fcro)]e 



(13) 



where si and ci denote the sine- and cosine-integrals, respec- 
tively. The model of a Plummer sphere has often been used 
in numerical simulations of the accretion and their eventual 
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disruption of satellite galaxies in massive parent galaxies. 
Plummer spheres have constant density cores, whereas nu- 
merical simulations of the formation of galactic haloes in 
cold dark matter cosm ology show that dark haloes ma y have 
a central density cusp (|Navarro, Frenk fc White! 19971 ) . Thus 
models of a Plummer or a Hernquist sphere should encom- 
pass the range of plausible models for satellite galaxies. The 
density in both models falls off radially steeper than found 
in the cold dark matter galaxy cosmogony simulations. This 
mimics the tidal truncation of satellite galaxies in the gravi- 
tational field of their parent galaxies. A more technical point 
is that a density profile as sh allow as r -3 as suggested by 
iNavarro. Frenk fc Whit el (1 19971 ) lends itself not easily to an 
analytic treatment, because the total mass would be diverg- 
ing without an outer cut-off. 



The Fourier transform of any potential with spherical sym- 
metry depends only on k = \k\. Thus it follows immediately 
from eq. (I19p that indeed the two acceleration components 



< v x > 



< v z > 







(20) 



as anticipated from symmetry reasons. Only in the direction 
of motion of the perturber there is a net effect. According to 
Newton's third law the drag force exerted on the perturber 
is given by rav = -nib < v y > where is the mass of a 
background particle, so that the drag force is anti-parallel 
to the velocity of the perturber. In order to evaluate the 
integrals over the wave numbers in eq. (|19p it is advanta- 
geous to switch from Cartesian form k x ,k y ,k z to a mixed 
representation k y , k = \/k x + k y + kf , arctan(k x /k z ), and 
we obtain for the deceleration the general result 



4 DYNAMICAL FRICTION 

The ensemble of stars is accelerated by the moving perturber 

as 



< v >= 



d 3 x / d 3 v/(x,v)V$i 



(14) 



where / denotes the full distribution function / = fo + /i . 
The contribution from fo cancels out, and introducing the 
Fourier transforms ([2]) we find 



< v >= 



d 3 x / d 3 v I d 3 kik$ k 



i[wt+k.x)] 
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(15) 



From symmetry reasons the acceleration vector < v > 
is expected to be oriented along the y-axis. In eq. (|15p 
the frequency uj, and similarly uj' , is given according to 
eqns. (01, ([TO} and (|13|) by ui — — k y vo — iX where we have in- 
troduced following Landau's rule a negative imaginary part, 
which we will let go to zero in the following. Moreover, 
$ k = $_ k so that the potential is a real quantity. Equa- 
tion (|15|) simplifies to 



d 3 v / d 3 kik$_ k / k e 2A ' 
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Using expression <(Sj this can be evaluated as 

ifcku 
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Next we observe that 
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so that 
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where erf denotes the usual error function. The Coulomb 
logarithm is defined as 



InA 



4tt 4 
G 2 m 2 



dk fc 3 |$ k | 2 



(22) 



In the case of a point mass formula J7J implies A — 
fcmax/fcm i n. Th is result was first obtained in this form by 
Kalnajs (1972) and is identical to Chandrasekhar's formula 
( Chandrasckhar 1943), if m + mi, fs m. The Coulomb loga- 
rithm diverges in the familiar way both on small and large 
scales, i.e. at fc~^ x and k~l , respectively. 

The Coulomb logarithm of the dynamical friction force 
exerted on a Plummer sphere can be calculated by inserting 
eq. dTOj into J22} leading to 



InA 



2 

: '0 



dk k K 1 (r\-,k) . 



(23) 



If the Plummer radius ro shrinks to zero, expression 
changes smoothly into the Coulomb logarithm of a point 



mass, because lim ro ^o roKi(kok) 



The integral over 



the square of the Bessel fun ctions in eq. (1231) can be eval u- 
ated using formula (5.54) of iGradshtevn fcRhvzikl (|2000D . 



2 i,2 

InA = -lofeia 



[A- 2 (r„fc m i n ) - Ko(rok mi n)K2(rok min )] (24) 



which is approximately 



InA 



- — ln(r fc n 



(25) 



in the limit of rofcmm «C 1. This modified Coul omb logarithm 
converges on small scales precisely as found bv lWhitg l| 19761 ) . 
but still diverges on large scales. A natural cut-off will be 
then the size of the stellar system under consideration. 

For a perturber with the density distribution of a Hern- 
quist profile we find a Coulomb logarithm of the form 



InA = / dk — [1 + rok cos(rofc)si(rofc) 
— rgk sin(rofc)ci(rofc)] 2 . 



(26) 



It can be shown using the asympt otic expansions of the 
sine- and cosine-integrals given by lAbramowitz fc Stegunl 
l|l972T l that the integrand in expression (|26l) falls off at large 
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Figure 1. Coulomb logarithms of the Plummer sphere (solid 
line), and a sphere with a Hernquist density profile (dashed line). 
The dotted line indicates — ln(27rro/A max )- ro denotes the radial 
scale lengths of the spheres and A ma x is the upper cut— off of 
the wavelength of the density perturbations (see text). The in- 
let shows the cumulative mass distributions of the Plummer and 
Hernquist models. 



k as 4ro(rofc) . Thus the Coulomb logarithm converges at 
small scales. This is expected because, although the density 
distribution has an inner density cusp, the deflecting mass 
'seen' by a field star with a small impact parameter scales 
with square of the impact parameter. At small wave num- 
bers a Taylor expansion shows that the square bracket in 
expression (|26[l approaches 1 so that we find a logarithmic 
divergence of the Coulomb logarithm as in the case of the 
Plummer sphere. 

In Fig. 1 we illustrate the Coulomb logarithms of the 
Plummer sphere and a sphere with a Hernquist density dis- 
tribution according to eqns. (|24[) and (|26|) as function of 
rofcmin- Fig. 1 shows clearly that at given mass small sized 
perturbers experience a stronger dynamical friction force 
than larger ones. Since the cut-off of the wave number at 
fc m in is determined by the radial extent of the stellar sys- 
tem, which corresponds roughly to one half of the largest 
subtended wave length A max = 27r/fc m i n , we use actually 
J"o/(Amax/4) as abscissa in Fig. 1. For comparison we have 
also drawn — in(rofemin) in Fig. 1. The insert shows the cu- 
mulative mass distributions of both mass models. The half 
mass radius of the Hernquist model measured in units of 
ro is about twice as that of the Plummer sphere. Thus for 
a proper comparison of the drag forces exerted on a Plum- 
mer sphere and a sphere with a Hernquist density profile 
the dashed line in Fig. 1 should be stretched by a factor 
of about 2 towards the right. But it is clear from Fig. 1 
that the drag force exerted on a Plummer is always larger 
that the drag on a sphere with a Hernquist density profile. 
This is to be expected because of its shallower density pro- 
file. The logarithm — £n(rofc m j n ), although being the asymp- 
totic expansion of the Coulomb logarithms (124[) and (|26[) for 
rofcmin — * 0, is not a good approximation at larger rofe m i n . 




Figure 2. Contour plot of m ■ InA for a Plummer sphere. The 
mass m is given in arbitrary units. The contour levels span the 
rnage from 0.5 to 4 at equidistant intervals of 0.5. 



There is a systematic off-set relative to the true Coulomb 
logarithms which is given explicitely in eq. (|25jl for the case 
of the Plummer sphere. 

So far we have treated the moving perturbers as rigid 
bodies. In reality, perturbers can be deformed by tidal fields. 
One source of the tidal field is the induced polarization cloud 
itself. However, its effect is expected to be small. The moving 
perturber induces the polarization cloud in the background 
medium which reacts back on the perturber. Thus the dy- 
namical friction force and any tidal fields exerted by the 
polarization cloud is of the order G 2 (cf. eq. 121) 1 . The de- 
formation of the perturber can be viewed as mass excesses 
and deficiencies. The momentum imparted by these 'extra' 
masses to the particles of the background medium will be 
of the order of G 3 and can be safely considered as a higher 
order effect. 

There is a further effect if the perturbers are grav- 
itationally bound systems themselves like globular clus- 
ters or dwarf satellite galaxies. Such objects can and do 
loose mass due to tidal shocking. This mass loss is pri- 
marily driven by tidal shocking due to the shrinking of 
the tidal radii at the inner pericenters of the orbits of 
the perturbing objects or by disc shocking when they 
pass through the galactic disc ([Gnedin fc Qstrikerl Il997l . 
iGnedin. Hernquist fc Ostrikerl Il999l ). Many studies, which 
cannot be all enumerated here, have shown that mass loss, 
which reduces dynamical friction, plays thus an important 
role for the evolution of the orbits of the satellite galaxies. 
In order to illustrate the effect of mass loss, on one hand, 
and the effect of a finite size of the perturber, on the other 
hand, we show in Fig. 2 the variation of the deceleration 
by dynamical friction exerted on a Plummer sphere as func- 
tion of the mass m and the radial scale length ro. As can 
be seen from the contour plot in Fig. 2 both effects can be 
of comparable magnitude. If one considers, for example, a 
perturber with a scale length of 0.01 (A ma x/4), a doubling 
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of its size has the same effect as a mass loss of 18 percent of 
the original mass. 

Finally we note that our analysis can be extended in a 
straightforwar d way, to anisotropic velocit y distributions of 
the field stars. iFuchs fc Athanassoula|[2005l have shown that, 
if the velocity ellipsoid is either prolate or oblate, the veloc- 
ity dispersion in the solution of the Boltzmann equation ((5} 
is replaced by an effective velocity dispersion which depends 
on the semi-axes of the velocity ellipsoid and its orientation 
relative to the wave vector k. This complicates the evalu- 
ation of the integrals with respect to the wave numbers in 
eq. (|19p considerably. We hope to address this problem in a 
forthcoming paper (Esquivel & Fuchs, in preparation). 
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